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ABSTRACT 

Type II supernovae (SNe), probably the most important contributors to stellar 
feedback in galaxy formation, explode within the very dense star-forming clouds, where 
the injected energy is most easily radiated away. The efficiency of type II SNe in 
injecting energy into the interstellar medium (ISM) and in re-heating a fraction of 
the original star-forming cloud is estimated with the aid of a two-phase model for the 
ISM of the cloud. We argue that when SNe explode the star-forming cloud has already 
been partially destroyed by ionizing light and winds from massive stars. SN remnants 
(SNRs) will first cause the collapse of most of the cloud gas into cold fragments, 
until the diffuse hot phase has a low enough density to make further radiative losses 
negligible. This is completed in ~3 Myr, with a modest energy loss of ~5 per cent of 
the total budget. We compute that a fraction ranging from 5 to 30 per cent of the 
cloud is reheated to a high temperature (from 10 5 to 10 7 K); these numbers are very 
uncertain, due to the very complicated nature of the problem, and to the uncertain 
role of thermal evaporation. Small star-forming clouds, less massive than ~10 4 M Q , 
will be destroyed by a single SN. In all cases, a high fraction of the energy from type 
II SNe (>80 per cent for large clouds, smaller but still significant for small clouds) 
will be available for heating the ISM. 

Key words: galaxies: formation - galaxies: ISM - ISM: bubbles - ISM: kinematics 
and dynamics 



1 INTRODUCTION 

The formation of stars and galaxies, as well as the state 
of the ISM, are regulated by the feedback processes related 
to the energy injection into the ISM itself by stars through 
winds, UV photons and SN explosions. In particular, a satis- 
factory model of galaxy formation requires that a significant 
fraction of the energy released by massive stars and SNe is 
given to the ISM and eventually to the hot, virialized gas 
component pervading the dark-matter halos. 

The greatest part of the stellar energy budget is pro- 
vided by type II SNe, associated to short-lived massive stars. 
The cosmological community often restricts to this feedback 
source alone, thus neglecting not only the contribution of 
UV and winds, that are produced by a subset of the stars 
that die as SNe and whose energy budget is less than the 
uncertainty in the energy of the single SNR, but also the 
contribution of type la SNe, associated to less massive stars. 
This way, energy is injected where young stars reside, i.e. in 
the star-forming regions, that are systematically the densest 
ISM regions in a galaxy. In this context the energy of SNe, 
that propagates into the ISM through blast waves, is very 
easily radiated away, giving rise to a low efficiency of energy 
injection. 



The computation of the energy lost by a SNR while it 
gets out of a star-forming cloud is not easy, due to a num- 
ber of complications. OB stars are highly clustered, both 
within the galaxy (they reside in the star-forming, molec- 
ular clouds) and within the star-forming cloud itself (they 
reside in associations). This way, most SNe explode in the 
hot bubbles created by previous explosions, giving rise to 
super-bubbles (SBs) more than isolated SNRs. Besides, the 
star-forming clouds are highly inhomogeneous, magnetized 
and dominated by supersonic turbulence. 

A further crucial point is that the first SNe explode 
from massive progenitor stars that have already pre-heated 
the surrounding ISM by UV light and winds. This process 
has been addressed by many authors, as for instance McKee, 
van Buren & Lazareff (1984), McKee (1989), Franco, Shore 
& Tenorio-Tagle (1994), Williams & McKee (1997), Matzner 
(2002), Tan & McKee (2004). Despite all of the uncertain- 
ties, some consensus is emerging on the fact that UV light 
and winds from massive stars are able both to sustain the 
observed level of turbulence in the star-forming clouds and 
to destroy them in a time of order of lOMyr. As a result, 
the bulk of SNe explode when the cloud is already in an 
advanced state of destruction. 

In a recent paper (Monaco 2004, hereafter paper I) we 
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proposed a model for feedback in galaxy formation in pres- 
ence of a multi-phase medium. In that model, the most dif- 
ficult piece of astrophysics to address was the efficiency of 
feedback in re-heating a fraction of the star-forming cloud 
to some high temperature, together with the amount of en- 
ergy lost by SNe before completing the destruction of the 
star-forming cloud. The first two quantities (fraction and 
temperature of re-heated matter) were left as free parame- 
ters, the last one (energy loss) was assumed to be negligible. 

In this paper we compute these quantities through a 
simple model for the explosion of SNe in a pre-heated star- 
forming cloud. The model assumes that SNe explode inside a 
cloud composed by a two-phase medium in pressure equilib- 
rium, with a hot diffuse phase and a cold phase fragmented 
into clouds. Due to the complexity of the problem, we are 
well aware that the model presented here is too naive either 
to represent the full complexity of the process or to give ac- 
curate predictions for the quantities involved. However, the 
results give insight on the physical processes in play and are 
informative enough at a qualitative level to constrain the 
order of magnitude of the quantities cited above; they can 
be used as a guide for future numerical simulations of the 
destruction of star-forming clouds. 

The paper is organized as follows. Section 2 presents 
the assumed initial conditions of the cloud when SNe start 
to explode, the model for the two-phase ISM, the fate of 
SNRs propagating in the cloud, the mass and energy flows 
within the components and the system of equations used. 
Section 3 presents the resulting predictions of the model. 
Section 4 discusses the case of small and dense star-forming 
clouds, where only a few SNe per cloud explode, and Section 
5 gives the conclusions. 

2 THE MODEL 

2.1 The state of the cloud at the first SN 
explosions 

Molecular clouds are dominated by supersonic turbulence 
(see, e.g., Solomon et al. 1987). According to recent simu- 
lations (Mac Low et al. 1998; Ostriker, Gammie & Stone 
1999; Mac Low 2003) supersonic turbulence in a compress- 
ible fluid decays over a few crossing times. This is true also 
in the case of magneto-hydrodynamical (MHD) turbulence. 
Besides, star formation takes place on longer times-scales 
and with a low efficiency, so that to obtain in a turbulent 
cloud a significant fraction of mass in stars (i.e. efficiency of 
star formation), comparable to the observed value ranging 
from 1 to 10 per cent (see, e.g., Carpenter 2000), it is neces- 
sary to sustain turbulence. Moreover, observations suggest 
that star-formation should not last more than ~10 Myr (see, 
e.g., Elmegreen 2002). 

UV light and winds from massive stars are likely re- 
sponsible for the destruction of the star-forming clouds. For 
instance, Matzner (2002) argued that the expanding HII re- 
gions are the most likely drivers of turbulence, and that the 
gradual photo-dissociation of H2 and the expulsion of re- 
heated material leads also to the destruction of the molecu- 
lar cloud, self-limiting the efficiency of star formation to the 
observed value. Requiring an equilibrium between the tur- 
bulence driven by the expanding HII regions and that dis- 
sipated by the turbulent cascade, and taking into account 




Figure 1. Mass flows in the model. 



the rate at which "blister" HII regions heat the gas to a 
temperature in excess of 10 4 K and expel it, he computed 
that a cloud will be destroyed in a time ranging from 10 to 
30 Myr, the lower values being valid for the largest molecu- 
lar clouds. The efficiency of star formation resulted ~10 per 
cent, nearly independent of cloud mass. His arguments are 
strictly valid for clouds with dynamical times shorter than 
the mean ionizing lifetime of massive stars (for a Milky Way 
ISM this amounts to M > 10 M@) and with escape veloc- 
ities smaller then twice the ionized sound speed (> 10 km 

According to this author, SNe do not contribute signif- 
icantly to the destruction of the cloud. In fact, the first SNe 
explode ~3 Myr after the formation of their progenitor stars 
and are indeed associated to the same OB stars that produce 
the HII regions. They are so massive that their remnants can 
go directly to the snowplow stage (see below) while still in 
the early free expansion stage; they end up being confined 
within the HII region, after having radiated away most of 
their energy. Their contribution to the momentum of the 
HII region is thus small. However, this argument is correct 
only for the first SN that explodes in each OB association. 
A second explosion within the same association would prop- 
agate into a rarefied hot bubble; in this case the blast would 
reach the ionization front before losing much energy. More- 
over, many SNe come from smaller stars, that explode later 
and are not necessarily associated to big HII regions. As a 
consequence, only the energy of the very first SNe will be 
completely lost; besides, energy injection by multiple SNe 
will become important later than 3 Myr, when the process 
of cloud destruction is already in an advanced state. 

While the precise origin of turbulence in star-forming 
clouds is still to be demonstrated, the energetic input of OB 
stars is likely to have a dramatic role in the destruction of 
the molecular, star-forming clouds. We thus deem it realis- 
tic to assume that when the first SNe explode the cloud is 
composed of two phases, a hot diffuse one, heated up by HII 
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regions, and a cold collapsed one, fragmented into cloudlets 
with a given mass spectrum. Pressure equilibrium between 
the two phases is assumed; this is justified by the finding of 
rough isobaricity in simulations of turbulent ISM (Kritsuk 
& Norman 2002; see Mac Low 2002 and Vazquez- Semadeni 
2002 for reviews). The expanding SNRs act in shaping the 
ISM within the cloud as follows (Figure 1): blasts in the 
adiabatic stage heat the diffuse phase, while in the snow- 
plow stage (see below) they collapse it. Thermo-evaporation 
of cold clouds within the expanding blasts transfers mass 
from the collapsed to the diffuse phase. Radiative cooling 
transfers mass from the diffuse to the collapsed phase. 

The initial conditions of the cloud are specified through 
its mass M c ioud and initial radius i? c ioud- Paper I shows that 
these quantities are related to the state of the ISM out- 
side the cloud; however, to keep the formalism simple we 
avoid making this connection explicit. In any case we assume 
that the external ISM is two-phase as well, with densities of 
"cold" and "hot" phases n c and n a respectively. Moreover, 
we assume for simplicity that the cloud is spherical. Follow- 
ing the discussion given above, a significant fraction of mass 
is put into the diffuse phase; the temperature of the collapsed 
phase is kept fixed to 100 K as in paper I. A fraction /* of 
the cloud is assumed to be in stars at the beginning of the 
calculation, which coincides with the instant at which the 
first SN explodes. Any further star formation is neglected. 

2.2 The two-phase medium within the cloud 

The state of the ISM, the expansion of SNRs and the mass 
and energy flows are modeled in a similar though sim- 
pler way as in paper I. Let Mdif and M co i be the mass 
in the diffuse and collapsed phases (with Mdif + M co i = 
(1 - /*)M c ioud), pdif and p co i their average densities (p = 
3M/47ri?j? loud ), T dif and T co i(=100 K) their temperatures, 
/idif and /icoi their mean molecular weights and fda and 
/coi their filling factors (with fdu + /col = 1). Moreover, let 
Fdu = Mdif/(M co i + Mdif) be the fraction of gas in the dif- 
fuse phase. The assumption of pressure equilibrium implies 
Pth/k = UdifTdif = n co iT co i, from which it is easy to obtain: 



Evaporative stage (t < t cv and t < t p ds) 



/col — 



1 + 



(1) 



-Fait Mdif T col 

(see Equation 2 of paper I). 

The collapsed phase is assumed to be fragmented into 
cloudlets, whose mass spectrum -/Vf rag is: 



A / "frag(mf r ag)dm fra g = iYo(mf ra g/l M©) draf ri 



(2) 



Here mf rag is the fragment mass in Mq and No a normaliza- 
tion constant such that p co i = J N{ Tag m{ rll gdm{ T!lg . At vari- 
ance with paper I we fix the exponent of the mass function 
to —2, which is a natural value in the presence of turbu- 
lence (see, e.g., Elmegreen 2001). The mass function is trun- 
cated below to a value mi set to 0.1 Mq (see paper I and 
the discussions below) and above by requiring unit prob- 
ability for the existence of at least one fragment (roughly 
m u = Mdoud/ ln(M c ioud/rai)). A typical radius <2f rag in pc 
is assigned to each cloud through the relation: 



fifrag = 0.104 /icolWcolOfrag 

Mfi . 



(3) 



In the following we will assume that the fragments are spher- 
ical. 
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Adiabatic stage (t BV < t < t p d s ) 
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PDS stage (t > t pdB ) 



R s (t) 


= -R pds (4t/3t pds - 1/3) 3 / 10 


V B (t) 


= v pds (4t/3t pds - 1/3) - 7 / 10 


Sth(t) t 


= E th (t pdB ) {0.398[1 - (t/1.169t') 14/5 ] 
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Skin(t) 


= Skin(t pd s)((^* P ds)Apd S ) _1/2 


M snp i(t) 


= M sw (t) [1- E th (t)/E th (t pdB )\ 



Table 1. Evolution of SNRs in the evaporative, adiabatic and 
PDS stages. R a , v a , f, M sw , M ev , i?th and -Bkin ar c respectively 
the radius, velocity, average temperature, swept mass, evaporated 
mass, thermal and kinetic energy of the SNR at a generic time 
t. The subscripts ev and pds refer respectively to the end of the 
evaporative stage and to the onset of the PDS stage. Here tg is 
time in units of 10 6 yr and Tq is Tdif in units of 10 6 K. All these 
quantities are valid for solar metallicity and the simple cooling 
function of Cioffi ct al. (1988) and paper I. t; (' = t pds exp(l), 
R' = R s (t') and the first term in parenthesis is present only for 
tpds < t < 3.16t pds (see Cioffi et al. 1988). 



Finally, we call v tur b the rms kinetic velocity of the ISM, 
eturb = 2u 2 urb the kinetic energy per unit mass and Pki n the 
resulting kinetic pressure. 



2.3 The fate of SNRs 

One star with mass > 8 Mq is formed each M* jSn Mq of 
stars (taken to be 120 Mq), so the total number of SNe is 
iV sn = /*M c ioud/M* iSn . The rate of SN explosions, R sn , is 
then computed as: 



Rsn 



/* Mdoud 
M* ian ilif e ' 



(4) 
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where tm c is the difference between the lifetimes of an 8 Mq 
star and the most massive star, assumed to be tm e =27 Myr. 
Any time-dependence of R sn is neglected. Each SNR injects 
-E51 10 51 erg of energy into the ISM. E51 is conservatively 
taken to be unity in the following, however its value is very 
uncertain and could be significantly higher. We call E sn (t) 
the total energy injected by SNe at the time t. We assume 
that SNe are homogeneously distributed within the cloud; 
this influences the estimate of the porosity of SNRs. The 
effect of the spatial clustering of OB stars will be commented 
later. 

Following McKee & Ostriker (1977) and paper I, the 
blasts are assumed to expand into the more pervasive dif- 
fuse phase; the cold dense fragments pierce the blast, that 
reforms soon after the passage. We will ignore the effect of 
the collapsed phase on the blast. At variance with paper I, 
we take into account the thermo-evaporation of cold clouds 
inside the expanding blasts, following the approach of Mc- 
Kee & Ostriker (1977) and Ostriker & McKee (1988); the 
saturation of thermo-evaporation is not taken into account. 
Assuming that the SNRs are not heavily mass loaded and do 
not lose thermal energy in the early free-expansion stage, the 
expansion of the SNRs is divided into three stages. For sake 
of clarity, we report the main properties of the evolution of 
the SNRs in Table 1. At the beginning thermo-evaporation, 
which depends on the 5/2-nd power of the average internal 
temperature, is very efficient, so that the interior hot gas 
is dominated by the evaporated mass. In this case the evo- 
lution of the remnant is given by the evaporative solution 
shown in Table 1, where: 

= . (5) 

In this equation the parameter <f> relates the actual thermal 
conduction to the Spitzer value (used to compute the rate of 
thermal evaporation), and its value depends sensitively on 
how effectively the magnetic fields quench thermal conduc- 
tion. Besides, spherical fragments present, at fixed mass and 
density, the smallest contact surface between the two phases, 
so any non-sphericity of the cloud will increase the effect of 
thermal evaporation; this can be mimicked by an increase 
of cp. Moreover, evaporation acts on the low-mass tail of the 
fragment mass function, and thus depends on the value of 
mi; keeping mi fixed, the uncertainty can be absorbed into 
4>. As a result this parameter is highly uncertain, and every 
value below or around unity is equally likely. The parame- 
ter q relates the blast speed to the average sound speed of 
the interior, and is in this case a 2 = 8. Finally, the quan- 
tity 1/al-ag is averaged over the mass distribution of clouds 
(Equation It is worth noticing that E has the dimension 
of a surface (pc 2 ), and that it diverges (E _1 vanishes) if 
thermal conduction is quenched. 

At later times the evaporated mass M cv becomes 
smaller than the swept mass M aw . From this moment we 
use the standard adiabatic solution, given again in Table 
1. Eventually, the interior gas cools and collapses into a 
thin cold shell that acts as a snowplow on the ISM. For 
the evolution in the Pressure-Driven Snowplow (PDS) stage 
we use the analytical model proposed by Cioffi, McKee & 
Bertschinger (1988), given again in Table 1, that fits rea- 
sonably well their detailed ID hydrodynamical simulations 
(confirmed by Thornton et al. 1998). 



The PDS stage can be reached when the blast is still in 
the evaporative regime. In this case we pass directly to the 
PDS solution of Cioffi et al. (1988), as the drop in the den- 
sity of the internal hot gas is very likely to quench thermal 
evaporation. 

The thermal energy of the hot interior gas in the evap- 
orative and adiabatic stages is respectively 55 and 72 per 
cent, the rest being kinetic. In the PDS stage the evolution 
of the thermal energy is given by Equation 3.15 of Cioffi et 
al. (1988) (reported in Table 1), while the kinetic energy is 
assumed to decay as ((t — t p ds)/tpds)~ 1,/2 ; this follows from 
assuming that all the mass and kinetic energy is in the ex- 
panding shell, so that i?kj n = M BVi v 2 /2, with M aw esc R^, 
v s oc Rs/t and _R S tx t ' 3 ; this is valid for t 2> t P ds- The 
interior mass is assumed to collapse into the snowplow at 
the same rate at which thermal energy is lost. 

The SNRs stop expanding in the following cases: (i) 
remnants stall by pressure confinement; this happens when 
the velocity of the blast equates the largest between the 
kinetic and the thermal velocity of the ISM. (ii) The poros- 
ity of the blasts, defined as Q = R Bn J Q Rl(t)dt/ R^ loud 
(and computed considering all the evolutionary stages of 
the SNRs), reaches unit values. The stopping time t s to P of 
the blast thus corresponds to the earlier of these events. 

The porosity given above is valid for a homogeneous 
distribution of SNRs. The non-uniform spatial distribution 
of OB stars will influence the estimate of Q if the radius 
-Rstop of the SNRs is small or comparable to the typical 
distance between associations. This is not the case; we have 
verified that in most cases the stopping radius of the SNR 
-Rstop is similar to the initial size of the cloud. This confirms 
the validity of the uniform distribution of SNRs as a first 
approximation. 



2.4 Mass and energy budget 

The onset of a PDS is important not only for its effect on 
the evolution of the blast but also for its effect on the diffuse 
phase, which is shocked to high temperature in the evapo- 
rative or adiabatic stage but compressed and cooled in the 
PDS stage. 

The rate at which the diffuse phase is swept by the ISM 
is M aw = Rs n M sw (tstop). Each SNR causes the following 
mass flows: (i) some collapsed mass M ev is evaporated to 
the diffuse phase; (ii) some diffuse mass M snp i is collapsed 
if PDS is reached. Both quantities are given in Table 1 and 
are always computed at the time tstop- In particular, the 
evaporated mass is given by the time-dependent term of Ta- 
ble 1 for t s top <t cv , while it remains constant after t cv (or 
tpds whenever it is smaller). The evaporation and PDS mass 
flows are then: 

Mov = R sn M cv (t stop ) , (6) 

M snp l = _R sn M snp l (tstop) • (7) 

Radiative cooling of the diffuse phase leads to a global 
decrease of thermal energy; in this process the density peaks 
cool dramatically and thus move to the collapsed phase, giv- 
ing rise to a cooling mass flow: 

M coo l = /cool ■ (8) 

£cool 
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As noted in paper I, the fraction / coo i depends on the de- 
tailed density structure of the diffuse phase, which is very 
difficult to predict without full-blown MHD simulation; no- 
tably, thermal conduction also influences it. f coo \ is thus left 
as a free parameter. The cooling time is: 



^cool 



3fcTdif 



n dif A(T di f 



(9) 



and is computed using the cooling function A(T) of Suther- 
land & Dopita (1993) 1 ; solar metallicity will be used 
throughout the paper. To avoid overcooling at the beginning 
of the integration a heating source is assumed to be present, 
such as to balance cooling for the diffuse phase present at 
the initial time. This is justified by the presence of the same 
UV photons responsible for the destruction of the cloud, and 
has an effect only during the first stages of evolution. 

The thermal energy of the diffuse phase -Edif is lost by 
radiation at a rate: 



E C ool 



Edit 



The diffuse phase gains energy by blasts at a rate: 

Eft, = RsnEth(tstop) , 



(10) 



(11) 



This includes the energy both in the evaporated and the 
heated gas. The rate of energy loss by snowplows is: 

3ft 



E 



snpl 



MmplTdif 



(12) 



where k is the Boltzmann constant and m p the proton mass. 
The kinetic energy per unit mass e tur b, measured in 



km^ 



of the collapsed phase decays through turbulent 



cascade as suggested by Mac Low (2002, 2003): 



e dis = 6 x f 7 u t 3 u 



r -1 , 2-2 

,L d km s yr 



(13) 



Here Ld is the driving scale of the turbulence (in pc) , taken 
to be twice the final diameter of the SNR (Mac Low 2002; 
Matzner 2002), or 4i? st0 p. The kinetic energy is continually 
replenished, provided that the coherent velocity field of the 
expanding SNRs is randomized by overlapping blasts. The 
fraction of kinetic energy that goes out of the cloud as a 
coherent velocity field is estimated as the ratio between the 
initial external area of the cloud (the one that contains all 
the stars) and the sum of all the areas of active blasts: 



T}2 

-"■cloud, 



Er s (i)dt. 



(14) 



As for the porosity, the integral is computed taking into 
account the evolution of R 3 (t) in all the stages up to t s t op - 
Of course / ex it is not allowed to exceed unity. The kinetic 
energy used to drive turbulence is then: 



en, = 5.03 x f0 7 (f - /exit)- 



-Rsn-Ekin 



km 2 s 2 yr 



(15) 



Mdif + Mcoi 

Here 5.03 x f0 r = f0 51 erg/(10 5 cm) 2 M Q . 

The thermal energy of the diffuse phase and the kinetic 
energy not used to drive turbulence are available for driving 
a SB into the external two-phase ISM with densities of "hot" 

1 We warn the reader that the analytical quantities given in Table 
I are computed using the very simple cooling function suggested 
by Ciofii et al. (1988) and used in paper I. 



and "cold" phases n n and n c (these phases have the same 
mean molecular weights as the diffuse and collapses ones). 
The exact modeling of the SB is beyond the interest of this 
paper, but the radius -Rdoud of the destroyed cloud will ex- 
pand with the SB. To follow this expansion we use some 
results that are valid in the adiabatic SB solution (Weaver 
et al. 1977; paper I). The velocity of the SB is assumed to 
be: 



v Bh (t) = ; 



L 



1/3 



R Bb 2/3 km s 



(16) 



where R s b is in pc and the mechanical luminosity L38, in 
units of 10 38 erg s _1 , is: 



Las = 3.16 x 10 5 (£ t h + /cxit^khO-Rs 



(17) 



The total energy used to drive the SB will be _E s b = J L38 dt. 
Moreover, we identify the radius of the shocked wind as the 
time- dependent radius of the cloud, -Rdoud, and set it to 
0.86ii s b. This way -Rdoud expands at a velocity: 



Vexp(t) = 0.86 « sb (i?cloud/0.86) . 



(18) 



The adiabatic expansion of the cloud leads to a further 
energy loss term for the diffuse phase. According to Weaver 
et al. (1977), a fraction 6/11 of the injected energy is given 
to the shocked external ISM. The energy lost by adiabatic 
expansion is then simply set as: 



E n 



Efb. 



(19) 



We warn the reader that the predicted Tdif depends sensi- 
tively on how the _E a d term is modeled. 

From all the mass and energy flows listed above the 
following system of equations can be written: 



Mdif = -M coo i - M snpl + M cv 

Mcol = Mcool + M snp l — Msv 

Edit = — -Bcool — E Bnp l + Efb — -Bad 

eturb = —edis + Cfb 



(20) 



This set of equations, together with Equation I18l for -Rdoud, 
is integrated with a standard Runge-Kutta code (Press et 
al. 1992). 



3 RESULTS 

The dynamical evolution of the system shows qualitative 
trends that are present for a very wide range of initial con- 
ditions and parameters. These trends are well illustrated by 
the example shown in Figure 2, relative to a cloud with 
M c i ou d=10 6 Mq, with initial conditions iW"dif=M co i=(l — 
/*)M c ioud/2, T dif =10 4 K and Wt U rb=10 km s 1 . Parameter 
values are set to E$i=l, <j)=0 and / C ooi=0.1. The initial ra- 
dius of the cloud is set using the same Equation |3] with 
the external "cold" phase density n c in place of n co i; for 
the external ISM we use as in paper I n c =10 and n n =I0 -3 
cm -3 , typical values for a galaxy disc. The figure shows the 
evolution of the masses of the two phases (Mdif and M co i), 
the energy released by SNe (E sn ) versus the energy used to 
drive the SB (E Bh ), the mass flows (M sw , M coo i, M snp i and 
M ev ), the thermal energy flows (-E coo i, E mp i, E th and E ad ), 
the state of the ISM (T di f, naif, n co i, Pth, -Pkin and v tU rb) 
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Figure 2. Evolution of a cloud with M c i out j=10 6 Mq for the reference choice of parameters given in the text, (a) Mass in the two phases, 
(b) Energy used to drive the SB compared to that injected by SNe. (c) Mass flows, including the rate at which mass is swept by SNRs. 
(d) Energy flows, (e) State of the ISM. (f) Mechanical luminosity and expansion of the SB. 
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Figure 3. The same as in Figure 2 for the example with <f> = 0.3. 
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and the quantities characterizing the expansion of the cloud 
within the SB (7? c i ou d, T38, Wox P and /exit)- 

At the starting time the cloud is dense enough to al- 
low SNRs to go into the PDS stage before being stopped. 
This creates a strong mass flow from the diffuse to the col- 
lapsed phase, with a dramatic drop of thermal pressure and 
densities of both phases. This increase of the mass of the 
collapsed phase does not imply a re-formation of the molec- 
ular cloud, as the fragments will be spread out within the 
volume of the expanding SB; they will later mix with the 
external cold phase. In this early stage pressure is domi- 
nated by the kinetic contribution due to the turbulent mo- 
tion of the cold phase. The porosity of SNRs soon reaches 
unity values, that are maintained throughout the evolution. 
When the diffuse phase is significantly depleted the SNRs 
merge before getting into the PDS stage; this happens in this 
example after 1.3 Myr. At this point energy is efficiently in- 
jected into the diffuse phase, whose temperature starts to 
increase. Once mechanical heating overtakes the effective 
heating term introduced in the equation, cooling becomes 
the dominant mass flow. This induces a further drop in the 
mass of the diffuse phase, which stabilizes after ~3 Myr to a 
constant value. From that moment the fraction of diffuse to 
collapsed mass and the density of the collapsed phase remain 
constant, while the other quantities follow the expansion of 
the cloud. Notably, the temperature of the diffuse phase in- 
creases gradually as a result of the efficient energy pumping. 
The energy used to drive the SB is significant after lMyr, 
energy losses are restricted to the first 3 Myr of evolution. 
Finally, the turbulent velocity iwb remains between 10 and 
25 km s~ 1 for the whole evolution. 

The introduction of thermal conduction changes the de- 
tails of the evolution but not its main properties. Figure 3 
shows the case with <j> — 0.3; the results depend only weakly 
on the value of 4> as long as it is comparable to unity. The 
snowplow and cooling mass flows are contrasted by evapo- 
ration, which however becomes dominant only after a few 
Myr. The process of collapse of the diffuse phase is com- 
pleted in some 3.5 Myr, but after that time the fraction of 
diffuse mass increases steadily. Due to the higher density of 
the hot phase, SNRs go to the PDS stage for the first 5.4 
Myr, with a consequent increase of energy losses. The struc- 
ture of the ISM is similar to before, but, due to its higher 
density, the diffuse phase is considerably colder and thermal 
pressure consequently lower. Due to the higher ndif, SNRs 
are confined at a smaller radius and / ex it never reaches unit 
values; this implies that some fraction of the kinetic energy 
is dissipated locally by turbulence. The evolution of the re- 
sulting SB is very similar, due to the weak dependence of its 
expansion velocity on mechanical luminosity (Equation ll6H . 

The system has been evolved with many combinations 
of initial conditions and parameter values. Figure 4 presents, 
as a function of cloud mass, the fraction of diffuse ISM Fda 
at 5 Myr, the temperature Tut of the diffuse phase at the 
same time, the time idestr required by the SNe to complete 
the destruction of the cloud (defined as the time at which the 
mass of the diffuse phase gets smaller than 5 per cent more 
than its minimum) and the fraction /i osrt of the total energy 
budget of SNe lost at 5 Myr. The quantities Tdif, Tdif and 
/i os t are measured at 5 Myr as at this time the destruction 
of the cloud is concluded in most cases; at later times Fdif 
and /iost are constant in absence of photo-evaporation, but 



increases otherwise, while Tdif increases. However, it must 
be noticed that the size of the cloud at this point is larger 
than the typical size of galaxy discs, so the applicability of 
this model at late times is doubtful. 

For the standard choice of parameters and initial con- 
ditions (as that used in Figure 2), from 3 to 18 per cent 
of the mass ends up in the diffuse phase, with Tdif ranging 
from 4 x 10 5 to 2.5 x 10 6 K; smaller clouds have lower R sn 
values, so they maintain more diffuse matter to a lower fi- 
nal temperature. These two trends tend to compensate each 
other, so the thermal energy contained in the diffuse phase 
is roughly proportional to the cloud mass. The time fdostr is 
~3 Myr and depends only weakly on the cloud mass, while 
the fraction of energy lost by radiation is well below 10 per 
cent but for the most massive clouds. 

When evaporation is switched on, the amount of diffuse 
matter increases significantly, especially for smaller cloud, 
while its temperature decreases to 10 5 K. This trend is 
stronger for larger 4> values, although the dependence on the 
actual <f> value is weak. As a result of the lower temperature, 
the thermal energy of the diffuse phase is lower than the no- 
evaporation case. Destruction times remain between 2 and 
5 Myr, while radiation losses increase but remain generally 
lower than 20 per cent. 

The other panels show what happens for the non- 
evaporative case when the parameter / coo i, the initial con- 
ditions (Fdit or Tdif) or the external densities are varied. In 
general, where more diffuse mass is obtained its temperature 
is lower, so the thermal energy in the diffuse phase does not 
vary as strongly as Tdif or Tdif. The destruction time in most 
cases ranges from 1 to 5 Myr, while the fraction of energy 
lost is quite stable. At variance with what found in paper I, 
the results of the integration depend very sensitively on the 
uncertain parameter / CO oi, and thus on the detailed density 
structure of the ISM. Moreover, they depend significantly 
on the state of the cloud when the SNe start to explode. Be- 
sides, the dependence on the external medium is not strong, 
though it is important to notice that with a denser external 
ISM Tdif tends to be lower and Tdif higher. 

A detailed justification of all the trends visible in Fig- 
ure 4 would shed some light on the dynamics of the system, 
at the cost of a lengthy discussion. The most important in- 
formation that this analysis yields is a direct impression on 
the robustness of the predictions. 



4 CLOUD DESTRUCTION IN THE 

ADIABATIC CONFINEMENT REGIME 

The system of Equations 1201 is valid as long as many SNe 
concur in determining the evolution of the system. More 
specifically, many SNe must explode in the first 3-5 Myr, 
which implies (for tm c =27 Myr) a total number of SNe ^>10 
and then a cloud mass M c i ou d^>10 4 Mq. For small clouds, 
the collapse of the diffuse phase and the onset of adiabatic 
confinement before PDS take place before the second SN 
manages to explode. In this case the present model simply 
does not apply. 

As shown in paper I, such small collapsing clouds are 
not found in spiral galaxies, where the Jeans mass (including 
non-sphericity and turbulent support) is rather high. This is 
a predicted result of the feedback regime in which SBs blow 
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Figure 4. i^difi ^dif (both at 5 Myr), tdestr an d /lost f° r many system as a function of cloud mass M c i ou( j. Standard refers to the choice 
of parameters of Figure 2, and is replicated in all panels. <j>: dotted, dashed and dot-dashed lines refer to </)=0.1, 0.3 and 1. f coo \. dotted 
and dashed lines refer to / coo i=0 and 0.3. Initial F^\f. dotted and dashed lines refer to F,jif=0.9 and F,jif=0.2. Initial T^if- dotted and 
dashed lines refer to 7dif=3 X 10 4 and 10 5 K. n c , nj,: dotted line refers to n c =100, 7ih=10 -2 cm -3 , dashed line refers to n c =l, ni!=10 — 4 
cm -3 . 



out of the disc, injecting their energy directly into the halo 
more than into the ISM. If the system is thicker or denser 
than a typical galaxy disc, SBs are kept pressure-confined by 
the hot phase, so they inject all of their energy into the ISM, 
which is then characterized by a much higher pressure and 
temperature of the hot phase. In this condition the Jeans 
mass is considerably lower, and collapsing clouds range from 
2000 to 10000 M©. As a consequence, only a bunch of SNe 
explode in each of them (this is especially true if /* is set to 
a low value). 

Let's consider the case of a M c i ou d=2000 Mq uniform 
cloud with density n c as large as 10 3 cm -3 , and let's call 
M2000 = Mcicmd/2000 M and niooo = n c /1000 cm" 3 , 
^cioud, computed from Equation |SJ will be 

-Rdoud = 2.5 M^n^oVj 3 pc . (21) 

Neglecting for the moment photo-evaporation and any pre- 



vious heating by UV photons, the SNR will go into the PDS 
stage well before exiting the cloud, as: 

R pds = 0.75 n^ooo 7 pc . (22) 

The blast will reach the outer boundary of the cloud after 

tout = 1.01 x 10 4 A^oJUioocf yr, (23) 

with a final velocity of 

v s = 69 M^Vooo 08 km s" 1 , (24) 

high enough not to be confined by the thermal or kinetic 
pressure of the cloud. At this point it will accelerate and 
fragment, spreading around the mass of the star-forming 
cloud. In this process it would lose nearly all of its energy. 
Any further SN would then explode in the rarefied bubble, 
reaching easily the outer boundaries of the cloud and merg- 
ing with the external hot phase. The fate of secondary SNe 
is very difficult to predict, but a rule of thumb would suggest 
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that while the energy of the first SN is lost, a significant frac- 
tion of the others will be injected into the medium. In this 
case the fraction of re-heated matter is likely to be rather 
low, in agreement with the extrapolation of the results of 
Figure 4. 

According to Matzner (2002), a single HII region can 
destroy a small molecular cloud. However, while UV photons 
can destroy all molecules and make the cloud unbound, thus 
lowering the average density of the cloud, it is likely that, for 
an external hot phase with P/k ~ 10 5 K cm -3 and Th ~ 10 7 
K, the region in which the first SN explodes will be much 
denser then the external hot phase. The SN will thus lose 
all of its energy and create a hot rarefied bubble in place of 
an overdensity 2 . So, the tentative conclusions given above 
remain valid also in the case of pre-heated clouds. 



5 DISCUSSION AND CONCLUSIONS 

Type II SNe are probably the most important source of stel- 
lar feedback, but they explode in the densest regions of the 
ISM, the star-forming molecular clouds. The energetic input 
by expanding HII regions and stellar winds can self-limit star 
formation to a low efficiency and destroy the cloud in ~10 
Myr. While the very first SNe will likely remain trapped 
within the HII regions, the bulk of them will explode when 
the cloud is already in the process of being destroyed. 

Under the assumption of a two-phase medium in pres- 
sure equilibrium, the system evolves into a configuration 
where most mass is in a collapsed, low filling-factor cold 
phase, while most volume is filled by a hot pervasive phase, 
able to confine the expanding SNRs while they are still in the 
adiabatic stage. This way most energy is efficiently pumped 
into the diffuse phase and used to power an SB expanding in 
the external medium. This is a very important point: thanks 
to the multi-phase nature of the ISM, only a few per cent of 
the energy injected by type II SNe is lost to radiation. 

The destruction of the cloud is completed by SNe in 
~3 Myr, with a loss of energy of ~5-10 per cent, to which 
one might add the energy of the very first SNe (roughly one 
per OB association) that are kept trapped within the HII 
regions. At the end of this process, a fraction of mass of the 
original cloud, ranging from 5 to 30 per cent, is heated to a 
temperature ranging from 10 5 to 10 7 K. This last quantity is 
especially uncertain; it depends sensitively on very uncertain 
parameters like cj>, f coo \ , on the specific initial conditions and 
also on the way energy is given to the external ISM to drive 
the SB. On the other hand, the total thermal energy given to 
the final diffuse phase is more robust. Anyway, the choices 
used in paper I of Fdif=0.1, Tdif=10 6 K and /i os t~0 are 
justified by these results. 

Excluding the very first explosions, SNe start to explode 
after a few Myr, and complete the destruction of the cloud 
in ~3 Myr. In this case, star formation cannot last more 
than several Myr. The destruction time inferred by Matzner 
(2002) of 10-30 Myr, based on the role of HII regions alone, 
can then be overestimated by a factor ~2-3. This would also 
apply to its predicted efficiency of star formation of ~10 per 

2 This can be checked by considering that i? p ds/^cloud ^ n iooo 
and v s oc n™^ 8 . 



cent. Low values of /*~5 per cent and of the duration of 
star formation, <10 Myr, are indeed in good agreement with 
observations (see, e.g, Carpenter 2000; Elmegreen 2000). 

These numbers refer to relatively large star-forming 
clouds, expected in galaxy discs where SBs are able to blow- 
out in the vertical direction. In this case, only 5-10 per cent 
of the energy budget of SNe is injected into the ISM (paper 
I), while a similar amount is lost in the destruction of the 
cloud. The rest of the budget, ~80 per cent, is injected into 
the halo, and is thus available to heat up the virialized halo 
gas, so as to preventing further cooling. 

Smaller and denser clouds, predicted in paper I to be 
associated to adiabatic confined SBs, are destroyed by one 
single SN. In this case the fraction of energy lost in the 
destruction is likely to be higher, while the fraction of re- 
heated mass is likely to be lower. 

Thermo-evaporation of cold clouds can play a signifi- 
cant role in this context. When <f> is not negligible the dif- 
fuse mass at the final time is more abundant and colder, and 
the energy loss is greater. However, thermo-evaporation does 
not change the qualitative behaviour of the system, and its 
contribution is still smaller than the uncertainties connected 
to the values of the other parameters. As a consequence, the 
main contribution to the reheating of the cold phase could 
be simply due to the residual effect of photo-evaporation by 
HII regions. 

While it is clear that only high-resolution MHD simu- 
lations of the ISM, able to include all the relevant physical 
processes in play, will be able to provide in the future robust 
quantitative predictions for the destruction of a star-forming 
cloud, the calculations presented here show that, thanks to 
the multi-phase nature of the ISM, a great part of the en- 
ergy from SNe will be able to leave the star-forming regions, 
thus being available to regulate the formation of galaxies. 
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